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Suzaku Observation of Strong Fluorescent Iron Line Emission from 
the Young Stellar Object V1647 Ori during Its New X-ray 

Outburst 

Kenji Hamaguchi^'^, Nicolas Grosso^'^, Joel H. Kastner^, David A. Weintraub^, Michael 

Richmond^ 

ABSTRACT 

The Suzaku X-ray satellite observed the young stellar object V1647 Ori on 
^ j 2008 October 8 during the new mass accretion outburst reported in August 2008. 

During the 87 ksec observation with a net exposure of 40 ks, VI 647 Ori showed 
a high level of X-ray emission with a gradual decrease in flux by a factor of 
5 and then displayed an abrupt flux increase by an order of magnitude. Such 
enhanced X-ray variability was also seen in XMM-Newton observations in 2004 
and 2005 during the 2003—2005 outburst, but has rarely been observed for other 
young stellar objects. The spectrum clearly displays emission from Helium-like 
iron, which is a signature of hot plasma {kT ~5 keV). It also shows a fluorescent 
iron Ka line with a remarkably large equivalent width of ~600 eV. Such a large 
equivalent width suggests that a part of the incident X-ray emission that irradi- 
ates the circumstellar material and/or the stellar surface is hidden from our line 
of sight. XMM-Newton spectra during the 2003—2005 outburst did not show a 
strong fluorescent iron Ka line, so that the structure of the circumstellar gas very 
close to the stellar core that absorbs and re-emits X-ray emission from the central 
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object may have changed in between 2005 and 2008. This phenomenon may be 
related to changes in the infrared morphology of McNeil's nebula between 2004 
and 2008. 

Subject headings: stars: formation — stars: individual (V1647 Ori) — stars: 
pre-main sequence — X-rays: stars 



1. Introduction 



Certain young stars dramatically increase their mass accretion rates by orders of mag- 
nitude (from 10~^Mq yr~^ to 10~^ Mq yr~^). These events are possibly triggered by thermal 
disk instabilities and are traced by dramatic increases in optical/IR luminosities; Such erup- 
tive pre-main sequence (PMS) stars are crudely classified as either FU Ori (FUor) or EX 
Lupi (EXor) types; the former are characterized by outbursts lasting decades, whereas the 
latter generally display sho rter outbursts, of duration a few months or years (see reviews in 
Hartmann &: Kenyorull996l ). Although PMS stars are known to be luminous X-ray sources, 
FUor or EXor events have been reported only rarely, such that the level of PMS X-ray 
activity during an outburst is poorly established. 

The young stellar object (YSO) V1647 O ri, which is deeply embedded in the L1630 
dark cloud {d ~400 pc, lAnthony-Twarod Il982l ) . underwent a strong optical/NIR outburst 
in 2003 December. This eruption afforded the first opportunity to record the sustained 
X-ray outburst of a rapid ly accreting YSO flKastner et al.ll2004l : iGrosso et al.ll2005l : iGrosso 



2006 



Kastner et al.ll2006l ). Multiple Chandra and XMM-Newton observations through this 



outburst demonstrated that the average X-ray flux varied in the same way and on approxi- 
mately the same timescale as the optical and IR brightnesses, although the X-ray flux also 
varied strongly (by a factor of up to 20) on timescales of less than a day. During outburst, 
the V1647 Ori X-ray spectrum had the characteristics of deeply embedded hot plasma (A^h 
~4.1xlO^^ cm~^, kT ~4.2 keV). This result appears to be explained best as star-disk mag- 
netic reconnection activity generated in association with the episode of very rapid mass infall 
jKastner et all l2004j . [2003 ). 



Following the V1647 Ori event, the EXor candidate VI 118 Ori recently displayed a 
mass accretion outburst; foUowup X-ray observations of this object detecte d a moderate 



enhancement in the X-ra y flux that was correlated w ith the optical/IR flux (lAudard et al. 



2OO5I : lAudard et aUboOoh . Iskinner et all J2006I . bood ) reported that two FUor systems with 
ongoing historical outbursts, FU Ori and V1735 Cyg, showed very hard spectra of kT'> 5 keV 
plasmas or equivalent despite relatively steady X-ray fluxes. Systematic surveys of low-mass 
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PMS (T Tauri) stars in th e Orion nebula and the Taurus dark cloud (e.g., iPreibisch et al. 
20051 : iTelleschi et al.l 120071 ) indicate that accreting T Tauri stars are less X-ray active (by a 
factor of ~2— 3 on average) than non-accreting T Tauri stars, suggesting that mass accretion 
activity may in fact suppress X-ray activity. Hence, the relationship between PMS mass 
accretion and high energy radiation remains unclear. 

V1647 O ri began a new optical/NIR ou t burst at a certain time between 2008 Jan. 2 
and Aug. 26 (lltagaki et al.l l2008t lAspinl l2008t lAspin et al.l l2009l ) . Our team later triggered 
an anticipated Target of Opportunity (ToO) observation of V1647 Ori with Chandra on 
September 18 for 20 ksec, in which the X-ray count rate was double the level observed by 
Chandra on 2004 March 7 f lKastner et al.ll2004l ) during the previous outburst in 2003—2005 
(Weintraub et al. in preparation). This new X-ray eruption offers us the first (and possibly 
unique) opportunity to measure X-ray emission during two outbursts from t he same YSO. 



We t herefore proposed a ToO observation with the Suzaku X-ray observatory (IMitsuda et al. 



20071 ). so as to obtain detailed X-ray spectral diagnostics and a ~day timescale light curve 
of V1647 Ori during its latest accretion outburst. With this information, we can derive the 
hottest plasma temperature and the intrinsic activity variation of V1647 Ori, unaffected by 
circumstellar absorption; we can also constrain models of the plasma geometry and of the 
response of neutral disk material to X-ray irradiation. This paper focuses on the results of 
the Suzaku observation. 



2. Observation 



Suzaku observed V1647 Ori for ~87 ksec on 2008 October 8, when the YSO's op- 
tical/NIR fluxes leveled off at a brightn ess similar to that observed in 2004 Feb.— Mar. 
fjOjha et al.ll2008l : lAspin fc Reipurthll2009l ). At the time of this observation, the observatory 
operated with t hree sets of telescope an d detector systems — the thin-foil X-Ray Telescope 



(XRT) modul e (jSerlemitsos et al 



Takahashi et al. 



2007 



20071 ) with the X-ray CCD detector (XIS: X-ray Imaging 



Spectrometer. iKovaraa et al.ll2007l) on its fo cal plane, and the Hard X-ray Detector (HXD: 



Kokubun et al.ll2007l ) with narrow collimator arrays. The XRT-I-XIS 



system focuses X-rays with a half-power diameter (HPD) of ~2' and covers the X-ray energy 
range up to ~10 keV. The XISl camera employs back-illuminated (BI) CCDs for the soft 
band sensitivity down to ~0.3 keV, while the other XIS cameras (XISO, XIS3) use front-side 
illuminated (FI) CCDs with sensitivity down to ~0.4 keV. The HXD detects hard X-rays 
above 15 keV. 



We started analysis from the cleaned event data in the distribution package, pre- 
processed with the version 2.2.11.22. The cleaned data excluded events recorded during 
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high background periods, mainly when the satelhte passed the South Atlantic Anomaly 
(SAA) 9 times during this observation, and earth occultations, each of which happened for 
~2.1 ksec in every Suzaku orbit (~5.4 ksec). The data did not suffer telemetry saturation. 
We thus selected events in good time intervals of all the XISO, 1 and 3 sensors. After this 
selection, the net exposure was 40,446 sec, i.e., ~46% of the observing time. For the XIS 
arf generation, we used the CALDB version released on 2009-02-03, to better reproduce the 
recent growth of contamination on the XISs, in particular XISO. The FI (XIS 0-1-3) and BI 
sensors collected 1216 and 624 counts between 1—8 keV, which include particle background 
and contamination from nearby sources by ~26% and ~39%, respectively. 

We generated time-averaged HXD spectra, following the standard analysis procedur^ 
and using the version 2.0 tuned backgrouncj^. The spectrum was consistent with the typical 
cosmic X-ray background spectrum and showed no hint of signal from V1647 Ori. 



3. Image Analysis 

Figure [1] shows an XIS color image of the V1647 field, combining all the XIS detectors 
(0-1-1-1-3). The XISs detected three bright peaks at the center, top, and bottom of the image 
and four faint sources at ~1.5' north-east, 4' south-east, 4' south, and 8' south-west from 
the central source. The bright sources to the north and south and the faint source to the 
south-east each have only one corresponding X-ray s ource in the 2002 C handra observation. 



respectively (source 14, source 8, and source 23 in ISimon et al.l |2004J ) within the Suzaku 



pointing uncertainty (<20")|f|. We measured their source positions from a 0.5—10 keV XIS 
image smoothed with the Gaussian function with cr=2 image pixel (16.7") by weighting 
photon counts within a 3x3 pixel region centered at the peak pixel. These p ositions were 



offset, on average, by (AR.A., ADec.) = (—8.3", —2.2") from the coordinates in lSimon et al. 



(120041 ). with a standard deviation of 2.0". We shifted the Suzaku image frame to match the 



Chandra coordinates. 

After these corrections, the absolute coordinates of the central bright source (a2ooo! (^2000) 
= (5''46'"13.35^ -00'^06'5.7") correspond to those of V1647 Ori within the 3a error range. 
Faint Chandra sources around V1647 Ori (9, 13, 15) are on the outskirts of the PSF and do 
not show any apparent local peaks at their locations in the Suzaku image. In the Chandra 



^http: / /heasarc. gsfc.nasa.gov/docs/suzaku/analysis/abc/ 
^ftp://lcgacy.gsfc.nasa.gov/suzaku/data/background/pinnxb_ver2.0_tuned/ 
'^suzakumcmo-2007-04 — http://www.astro.isas.ac.jp/suzaku/doc/ 
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observation in September 2008 (ObsID: 9915), these sources were far weaker than (had ^3% 
of the photon counts of) V1647 Ori. Their contamination of V1647 Ori in the Suzaku data 
would therefore appear to be minimal. Th e soft source to the north-east from V1647 Ori 
is probably a blend of sources 20 and 21 in ISimon et al.l (j2004f ). while the very hard source 
to the south is a blend of sources 7, 10, and 11, which are members of SSV 63, a cluster of 
class I protostars. 



4. Timing and Spectral Analysis 

The north-east source heavily contaminated the V1647 Ori region below ~1 keV, but its 
contamination was negligible above 4 keV. To minimize the contamination and maximize the 
photon statistics, we defined a source region with a 1.5' radius circle, excluding the circular 
region within 40" of the north-east source (see the right panel of Figure [1]) . To cancel out the 
remaining contamination, we defined a background region that is symmetric to the source 
region with respect to the north-east source. The background-subtracted light curves and 
spectra presented in Figure [SHU were generated from the events within these regions. 

Figure |2] shows the background subtracted X-ray light curve of V1647 Ori covering the 
energy range 1—8 keV, combining all the XIS (0+1+3) data, and using a bin size of 2 ks. 
The few missing data points are time periods where few good-time intervals (<1%) were 
available. The count rate varied strongly, by a factor of >18, during the observation (count 
rate of the maximum bin: ~0.11 counts s~^, la upper limit count rate of the minimum bin: 
~6xl0~'^ counts s~^). The X-ray light curve shows a gradual decrease, to a level consistent 
with zero flux, over the first ~60 ks of the observation, then a sharp increase in flux at 
~60 ksec followed by a decay with significant spikes and dips. 

We produced background subtracted soft band (1—4 keV) and hard band (4—8 keV) 
light curves and calculated the hardness ratio (HR) defined by {H — S)/{H + S), where S 
and H are the soft and hard band count rates, respectively (bottom panel of Figure The 
HR was ~ —0.5 between 10 ksec and 60 ksec into the observation and increased to ~0.0 
after ~60 ksec, remaining at that level for the rest of the observation duration. The abrupt 
flux increase at ~60 ksec, accompanying an HR increase, may be explained by an increase in 
kT (see also spectral analysis in the latter paragraphs), while fluctuations in flux after that, 
which are apparently uncorrelated with the HR, caused by changes in the plasma emission 
measure. Such flux and HR variations were seen in the XMM-Newton observation in 2004, 



as well ( iGrosso et al.ll2005l ) 



We merged spectra of the FI sensors and produced time-averaged FI and BI spectra 
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(Figure |3]). All the spectra in our analysis were binned so as to obtain at least 20 photon 
counts per bin. We ignored energy bins below 1.4 keV, where contamination from the north- 
east source exceeds the net signal from V1647 Ori. The resulting spectra reveal a strong 
fluorescent iron line at 6.4 keV from quasi- neutral iron atoms, as well as the helium-like iron 
line at 6.7 keV, which is a signature of hot plasma. These features can be seen clearly in 
the Fe Ka line region of the merged XIS spectra (right panel of Figure [3]). The spike at 
5.9 keV possibly originates from contamination of the ^^Fe calibration source scattering over 
the CCD chips. 

The spectrum is well fit by a model consisting of a single-component, optically thin- 
thermal plasma emission (APEC) combined with a narrow Gaussian line with the center 
energy fixed at 6.4 keV for the iron fluorescent emission, suffering photoelectric absorption 
by neutral gas (wabs). The best-fit parameters are listed in the model "i" in Table [TJ 
The best-fit plasma temperature (~4.1 keV), elemental abundance (~0.50 solar), hydrogen 
column density (~3.6xl0^^ cm"^), and the average observed flux (~4.4 xlO~^^ ergs cm~^ s~^ 
between 0.5—8 keV) are similar to tho se inferred for V1647 Ori during the XMM-Newton 



observation in 2004 (IGrosso et al.ll2005l ). In the 2.8—8.0 keV band, the absorption corrected 
X-ray luminosity during the Suzaku observation (log Lx ~30.9 ergs s~^) was twice that seen 
during the XMM-Newton observation (log Lx ~30.5 ergs s~^). One remarkable difference, 
however, is the firm detection of the iron fluorescent line at 6.4 keV. The equivalent width 
(EW) in 2008, ~600 eV, was about a factor of 6 higher than that found from the marginal 
detection of the iron fluorescent line during the XMM-Newton observation in 2004. 

Because the spectrum changed so dramatically near the midpoint of the Suzaku expo- 
sure, we divided the observation into two phases, before and after 56 ksec from the observa- 
tion start, as shown at the top of Figure [2l We simultaneously fit these spectra with their 
elemental abundances tied. The best-fit results are shown in the model "ii" in Table [T] and 
Figure m The plasma temperature was significantly lower in the early phase, while the other 
parameters, especially the hydrogen column density, do not appear to vary significantly. 
This would suggest that the plasma heated up and/or that a very hot plasma emerged at 
~56 ksec after the start of the observation. The best-fit EW of the iron fluorescent line 
was huge, close to 1.5 keV, during the "early" phase though the error was large enough to 
include the best-fit EW in the "late" phase. The model "iii" in Tabled] with their absorption 
and abundance parameters tied gives a similar result. Such large fluorescent iron line EWs 
have been observed from spatially extended X-ray reflection sources that do not contam inate 



emission from their irradiating sources (e.g., iKoyama et al.lll996t ICorcoran et al.l 120041 ). 
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Discussion 



The Suzaku observation was performed between ~40 and ~280 days after the onset of 
the new accretion outburst from V1647 Ori. The first (2004) XMM-Newton observation of 
V1647 Ori during its 2003-2005 outburst was performed with similar timing, ~4 months 
after outburst onset, with net exposure time (37—39 ksec) similar to that of the Suzaku ob- 
servation. Since these two observatories have similar effective areas in the 0.5—10 keV energy 
range, these observations are well suited to comparing the X-ray properties of V1647 Ori 
between the 2003 and 2008 outbursts. 

The flux of X-ray emission from V1647 Ori strongly varied, by at least a factor of 18, 
during the Suzaku observation. The flux bottom is not well constrained given the contam- 
ination from the north-east soft source, but the XMM-Newton observation in 2004 showed 
a similar range of variation of between 0.5—12x10"^^ ergs cm~^ s"^. Both observations 
showed similar abrupt flux increases with correspondingly HR increases and no significant 
HR variations after the increases. The plasma parameters derived from the Suzaku spec- 
tra were sim i lar to the best-fit result of the XMM-Newton spectrum (model 1 in Table [H 
Grosso et al.l 120051 ). These similarities suggest that the new outburst in 2008 was driven 
by a mechanism similar to that of the outburst which started in 2003 October and lasted 
for ~2 years. Similarities between the 200 3—2005 and latest outbursts are also seen at 
near- infrared wavelengths (lAspin et al.ll2009l ). 



While VI 647 Ori displayed strikingly similar X-ray activity patterns during its two mass 
accretion outbursts, there appears to be a very wide range of X-ray behavior and pre/post- 
outburst X-ray luminosities among the other FUor and EXor systems observed in X-rays 
at the onset of and/or during eruptions. The FUor systems FU Ori and V I 735 Ori emit 



strong X-ra ys at levels of log Lx ~31 ergs s ^ from hot plasma, kT > 5 keV ( ISkinner et al. 



2006 



20091 ). similar to that seen from V1647 Ori. On the other hand, several FUor and 



EXor systems did not show rem arkable X-ray activi ty associated with their mass accretion 
outbursts (EXor — VI 118 Ori: kudard et allbood FUor — V1057 Cyg and V1515 Cyg: 
Skinner et al.l 120091 ). These particular X-ray characteristics are apparently not directly re- 
lated to the optical outburst criteria that are used to classify objects as FUors or EXors (e.g., 
outburst duration). The X-ray behavior during outburst therefore might be more reasonably 
connected to the other YSO characteristics, such as rotational velocity and/or strength of 
the magnetic field. 

The Suzaku spectra in 2008 showed an impressively strong fluorescent iron Ka line 
feature at 6.4 keV. Such a l arge EW value has bee n repor ted only from a few young stel- 
lar obiects (X^ or IRS 7B: iHamaguchi et al.l 120051 IRSl: Iskinner et all 120071 V1486 Ori: 
Czesla fc SchmittI 120071 ). The observed iron fluorescent EWs are more than 10 times as 
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large as that for a 4 keV plasma ir radiating an infin ite plane of solar iron abundance in 
a Monte-Carlo simulation (^60 eV, iDrake et al.ll2008l ). so that special physical conditions 
would be required to explain such a large EW. Plausible conditions that can elevate the 
EWs are (1): supersolar iron abundance, (2): re-emission by a neutral absorber, e.g., a 
flaring disk and/or infalling envelope, (3): time-delay effect of fluorescence intensified by 
radiation from a magnetic flare, (4): line-of-sight obscuration of the irrad iating plasma, (5) : 
excitation by non-therm al electrons produced by a magnetic flare (e.g., lOsten et al.l 120071 : 
Czesla fc SchmittI 120071 ). Among them, the conditions (3) and (5) would require a big flare 
to produce the large EW, and which was not observed during the Suzaku observation. The 
condition (2) can only increase the EW by a factor of 2. The condition (1) is unlikely be- 
cause, in this condition, the photosphere and/or disk abundances would have had to increase 
by a factor of six between 2004 and 2008. We thus favor the condition (4) as the most likely 
explanation. 

Assuming the presence of a reflector with solar iron abundance on an infinite plane, 
~90% of the plasma has to be blocked from our direct view by an optically thick absorber 
in order to explain the large iron fluorescent line EW measured for V1647 Ori in 2008 by 
Suzaku. Such a geometry is plaus ible if most of the plasm a is hidden be hind the stell a r core 
and/or disk (see the discussion in lHamaguchi et al.ll2005[ ). We note that lAspin et al.l (|2009[ ) 
demonstrated a change in the r' band morphology of the reflection nebula associated with 
V1647 Ori (McNeil's Nebula) between 2004 and 2008, and suggested that dust obscuration 
close to V1647 Ori may play an important role in defining the observed morphology of the 
nebula. The circumstellar gas structure very close to the stellar core may have changed 
between 2004 and 2008, and this change may have altered both the reflection paths of X-ray 
emission and the openings in visible light between the star and McNeil's nebula. 



We greatly appreciate the Suzaku review committees positive response to our DDT 
request. This work is performed while K.H. was supported by the NASA's Astrobiology 
Institute (RTOP 344-53-51) to the Goddard Center for Astrobiology. J.K.'s research on 
X-rays from erupting YSOs is supported by NASA/GSFC XMM-Newton Guest Observer 
grant NNX09AC11G to RIT. 
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Table 1. Spectral Fits 



Model 


Phase 


kT 


Z 




Line Flux" 


Line EW'' 


log Ly," 


[d.o.f.) 






(keV) 


(solar) 


(10^2 cm-2) 


(10-6 ph cm-2 s-i) 


(eV) 


(ergs s^-"-) 




i 


whole 


4.1 (3.2-5.6) 


0.50 (0.31-0.73) 


3.6 (2.9-4.4) 


2.7 (1.6-3.8) 


588 


31.1/30.9/31.3 


1.03 (85) 


ii 


early 


1.9 (1.4-2.7) 


0.65 (0.36-1.02) 


5.1 (3.7-7.1) 


2.1 (0.66-3.6) 


1328 


31.3/30.7/31.4 


1.00 (84) 




late 


5.3 (3.8-7.5) 


=early 


3.6 (2.9-4.5) 


4.5 (2.3-6.8) 


549 


31.2/31.2/31.5 




iii 


early 


2.3 (1.8-3.2) 


0.62 (0.32-0.97) 


4.0 (3.3-4.9) 


2.4 (0.33-3.5) 


1355 


31.1/30.7/31.2 


1.02 (85) 




late 


4.7 (3.5-6.8) 


=early 


=early 


3.8 (2.1-7.0) 


459 


31.2/31.2/31.5 




1^ 


2004 April 4 


3.0 (2.4-3.9) 


0.8 (0.5-1.3) 


2.9 (2.5-3.4) 






30.8/30.5/31.0 


1.11 (129) 



Note. — "Total photon flux in the narrow Gaussian line with a fixed center at 6.4 keV. ''Measured with the xspecl2 command "cqw range C 
Gaussian^componen t^ID'^ . '^Absorption corrected luminosity between 0.5—2.8/2.8—8/0.5—8 keV, assuming d =400 pc. "^Model #1 in Table 1 oi 
ICrosso et al.l l|2005l V The parentheses between the second and fifth columns show 90% confidence ranges. 
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Fig. 1.— Left XIS 0+1+3 color image of the V1647 Ori field (red: 0.5-2 keV, green: 
2—4 keV, blue: 4— 10 keV) after correcting the absolute coordinates using a Chandra result 
( ISimon et al.ll2004j ). The image was smoothed by the Gaussian function with a =3 image 
pixels (~25") and displayed in a linear scale for the 5 counts pixel"^ range in each color, 
offset at the background level determined from a source free region. The extended blue 
areas at the 3 detector corners are ^^Fe calibration sources. Green circles with numbers and 
a black circle at the center (V1647 Ori) are X-ray sou r ce positions detected in the Chandra 
observation on 2002 November 14 fjSimon et al.l l2004l : iKastner et al.l 120041 ) . Radii of these 
circles are proportional to Chandra photon count rates between 0.5—10 keV in logarithmic 
scale, as shown at the bottom right corner. The grey dots show the Chandra ACIS-S FOV. 
Right: monochromatic XIS 0+1+3 unsmoothed image of the same field between 0.5—10 keV. 
The cyan and red circles show source and background regions, respectively; the region within 
the magenta circle was excluded from these source and background regions. 
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Fig. 2. — Top panel: Background subtracted light curve of V1647 Ori between 1—8 keV, 
produced from the XIS 0+1+3 data. Each bin has 2 ksec. Bottom panel: The hardness 
ratio curve with 4 ksec bins defined by {H — S)/{H + 5), where H and S are the hard 
(4—8 keV) and soft (1—4 keV) band count rates, respectively. The label above the top panel 
shows intervals defined for phase-resolved spectral analysis. The origin of the x-axis is 2008 
October 8, 14'^33"n0^ in UT. 
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Fig. 3. — Left Top: XIS FI (black) and BI (red) spectra of V1647 Ori. The solid lines 
show the best-fit result of the spectra by a model of 1-temperature thin-thermal plasma 
emission (APEC) with a Gaussian function at 6.4 keV for the iron fluorescent line, suffering 
photoelectric absorption by quasi-neutral gas. Left Bottom: residuals of the values from 
the best-fit model. Right: XIS 0+1+3 spectrum near the iron K line region. The vertical 
axis shows the ratio of the emission against the continuum in the best-fit model (model "i" 
in Tabled]). Important emission lines are marked by lines at their rest energies. 




Fig. 4. — X-ray spectra of the early (left) and late (right) phases (model "ii"). 



